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The phase-ordering dynamics that result from domain coarsening are considered for itinerant 
quantum ferromagnets. The fluctuation effects that invalidate the Hertz theory of the quantum 
phase transition also affect the phase ordering. For a quench into the ordered phase there appears a 
transient regime where the domain growth follows a different power law than in the classical case, and 
for asymptotically long times the prefactor of the t 1 ^ 2 growth law has an anomalous magnetization 
dependence. A quench to the quantum critical point results in a growth law that is not a power-law 
function of time. Both phenomenological scaling arguments and renormalization-group arguments 
are given to derive these results, and estimates of experimentally relevant length and time scales are 
presented. 

PACS numbers: 75.20.En; 75.40.Gb; 05.70.Ln; 73.43.Nq 



I. INTRODUCTION 

When a many-body system capable of a phase transi- 
tion from a disordered phase to a phase with long-range 
order is suddenly taken, by changing one or more param- 
eters, from the disordered phase to the ordered one, an 
interesting question is how the long-range order will de- 
velop as time goes by and the system approaches equilib- 
rium. Such a sudden transformation is called a "quench" , 
and the phase-ordering dynamics after the quench can be 
studied by kinetic methods similar to those used for the 
critical dynamics near the phase transition^ The quench- 
ing problem is of broad interest since it is applicable to 
a large variety of physical systems that undergo phase 
transitions, ranging from magnets to liquid helium to the 
early universe £ The phase ordering occurs by means of 
the growth of domains that arise from spontaneous fluc- 
tuations, and the linear size L of these domains^ obeys a 
power-law as a function of time t for sufficiently large t: 
L(i) oc t x l z , with z a dynamical exponent^ In addition, 
the pair correlation function is observed, both experimen- 
tally and numerically, to obey a simple scaling law 



principles has been given. This phenomenology has so 
far been applied to classical systems, but there is no rea- 
son to expect that it will not be valid for quantum sys- 
tems as well. Quantum phase transitions are known to 
differ in crucial aspects from classical ones*££ and one 
needs to ask whether these differences affect the phase 
ordering properties as well. In this paper we investigate 
this problem for the case of a quantum ferromagnet and 
show that the phase ordering kinetics are indeed affected 
in dramatic ways. 



II. REVIEW OF RESULTS FOR CLASSICAL 
MAGNETS 

In order to motivate our approach and put it into con- 
text, we first briefly recall the known results for phase 
ordering in classical magnets. The dynamical equation 
that governs the time evolution of the order parameter 
in an isotropic Heisenberg ferromagnet isl 



CM) = ( < p(x,t)-cl>(0,t)) = f(r/L(t)), 



(1.1) 
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Here r = \x\, / is a scaling function, <f> is the order 
parameter field, and (. . .} denotes a statistical average. 
We take the order parameter to be a real 3-vector, that 
is, we consider Heisenberg magnets. 

The facts stated above, although well established^ are 
purely phenomenological; so far no derivation from first 



Here A is a spin transport coefficient, and 7 is a gyromag- 
netic ratio. The Langevin force £ is of no consequence 
for the problem of phase ordering and we will neglect 
it.— H is the Hamiltonian or free energy functional that 
governs the equilibrium properties of the system. A clas- 
sical isotropic Heisenberg ferromagnet is described by a 
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/> 4 -theory with parameters r, c > 0, and u > 8 
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(2.2a) 



In Fourier space, this corresponds to a Gaussian vertex 

r(fc) = r + cfc 2 . (2.2b) 

The first term in Eq. (|2.1[) describes dissipative dynam- 
ics for a conserved order parameter; this is Model B in 
Ref.fl Including the second term takes into account the 
precession of spins in the effective magnetic field created 
by all other spins; this is Model J in Ref. [l|. 

The phase-ordering problem for Model B (7 = 0) has 
been studied by a variety of analytic techniques as well 
as by simulations A The result is a dynamical exponent 
z = 4; i.e., the linear domain size grows for long times as 
L(t) oc i 1 / 4 , with a prefactor that is independent of the 
equilibrium magnetization tuq. This result is plausible 
from a simple power-counting argument: Eq. (|2.1|) has 
the structure of a continuity equation for each component 
(j) a of the order parameter, <9 t <p a = — V-j'q., and for power 
counting purposes the current j a can be identified with 
the domain growth velocity dL/dt times itiq. Assuming 
that each gradient can be identified with a factor of 1/ L, 9 
this leads to dL/dt ~ cX/L 
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or 



L(t) oc (cA) 1/4 i 1/4 , (Model B). (2.3) 

The additional term in Model J (7 ^ 0) describes spin 
waves whose dispersion relation follows from Eq. (|2.1|) : 

u(k) = D(m )k 2 , (2.4) 

with D(mo) = 7 c mo = c£1l, where 51l = 7 mo is 
the Larmor frequency related to the equilibrium mag- 
netization m,Q. This is consistent with results obtained 
from microscopic models^ For the phase ordering prob- 
lem, Model J was studied in Ref. [l3|. The same power- 
counting arguments as for Model B above suggest 



dL 
~dt 



cX 



~L 



(2.5) 



According to Eq. (|2.5|) , the time dependence of L will 
cross over from the t 1 ^ 4 behavior characteristi c for M odel 
B to a t 1 / 2 behavior at a length scale L\ = \/ A/Ol, 



L(t) oc 



(cA) 1 / 4 ^/ 4 ifL<Li 



(Model J). 



(2.6) 

A numerical solution of the dynamical equation was 
found to be in good agreement with this expectation.— 
This concludes our review of known results for classical 
magnets. 



III. QUANTUM FERROMAGNETS 

A. Mode-mode coupling effects 

The above results hold if the equilibrium properties 
of the ferromagnet are described by Eqs. (|2.2[) , and if 




FIG. 1: Schematic phase diagram in a temperature (T) - 
pressure (p) plane with a phase separation line separating 
a paramagnetic (PM) phase from a ferromagnetic (FM) one. 
Shown are the quantum critical point (QCP), the classical 
critical regime (I), and static (IIa,b) and dynamic (III) quan- 
tum critical regimes. The arrows denote a quench into the 
ordered phase (A), and a critical quench (Bh respectively. 
See the text for further explanation, and Ref. [g for a general 
discussion of magnetic quantum criticality. 



there are no other soft modes that couple to the order 
parameter. A counterexample is phase separation in bi- 
nary fluids, where one needs to take into account that 
the local fluid velocity contributes to the order parame- 
ter transport^ The net result is equivalent to a nonlocal 
free energy, or dynamic equation, and this is obtained 
explicitly if the additional soft modes are integrated out. 
At low temperature (T) in itinerant fcrromagnets a sim- 
ilar phenomenon occurs; viz., a coupling of the order 
parameter to soft particle-hole excitationsj 14 i 15 These 
mode-mode coupling effects invalidate Hertz's mean-field 
theoryji£ and they change both the critical behavior at 
the ferromagnetic quantum phase transition^ and the 
magnetization dependence of the magnon dispersion rela- 
tion in the ordered phased Here we investigate how these 
effects influence the phase ordering following a quench 
into the quantum regime, which can be realized by a 
pressure quench at fixed low temperature in a system 
where the ferromagnetic quantum phase transition can 
be tuned by hydrostatic pressure, see Fig. [TJ Examples 
of such systems include UGe^ and MnSii2£ The equi- 
librium quantum phase transitions in these systems have 
been studied experimentally in some detail, and our pre- 
dictions for the phase ordering should be amenable to 
experimental checks using similar methods. 

The nature of the mode- mode coupling effects depends 
on whether the system is dirty or clean, i.e., whether or 
not quenched disorder is present. In the clean case, they 
lead to a fluctuation-induced first-order transition^ so 
the magnetization cannot be made arbitrarily small. In 
the dirty case, the quantum phase transition generically 



3 



is of second order. We will focus on the latter j22 where 
the size of the quantum effects we consider is not lim- 
ited by a nonzero minimum value of the magnetization. 
Quenched disorder also ensures that the transport coef- 
ficient A will remain finite even at T = 0. The form of 
the dynamical equation (|2.1|) will thus not be modified 
by quantum mechanics. We will briefly discuss clean sys- 
tems later. 

Effectively, the mode-mode coupling effects in a dirty 
system in 2 < d < 4 dimensions lead to a non-local free- 
energy functional that contains a V d_2 term in addition 
to the usual V 2 term, or to a Gaussian vertej*^^ 



r(fe) = r + £\k 



id— 2 



■k 2 



instead of Eq. (|2.2bj) . with c > 0. In the ordered phase, 
this nonanalyticity is cut off by the magnetization, with 
too ~ k 2 . It is also cut off by T > 0. We will discuss 
the latter effect, as well as the behavior at criticality, be- 
low; for now we assume any length scale associated with 
temperature to be the largest scale in the system, which 
effectively sets T = 0, and a quench well into the or- 
dered phase, see trajectory A in Fig. [TJ 2 ^ Let us denote 
the length scale where mo cuts off the nonanalyticity by 

— 1/2 

£2( m o) oc to , and the length scale beyond which 
the nonanalyticity dominates by L* — (c/c) 1 ^ 4- ^. (We 
will determine and discuss L 2 and L* in more detail be- 
low.) Since the dissipative term and the torque term 
in Eq. (|2.I[) are both proportional to 5H/5(f>, they are 
equally affected. Equation (|3.f [) implies that, effectively, 
5H/5(j> = —cV 2 cf> which would result from Eq. (|2.2a[) 
gets multiplied by a function /(V, mo), where V stands 
for the appropriate inverse length scale, which is |fe| for 
the magnon dispersion, and 1/L for the phase ordering 
problem. For scales larger than L* , we have 

.,_ , \{L* V) d - 4 cx V d - 4 iff/V<£ 2 
/(V, too) oc I {l Vi2)d _ 4 ^ m _ (4 _ d)/2 . f > ^ 

(3.2) 

For an illustration of these effects, let us consider the 
magnon dispersion relation. The above considerations 
result in Eq. (|2.4p with a modified _D(too), viz., D(mo) oc 
TO ( d -2)/ 2 ^ f or |^| ^ i/£ 2j and i n W (A;) oc |fc| d ~ 2 for \k\ ^> 
1/L 2 . The former result was first obtained in Ref. [TH 
from microscopic considerations, and we have reproduced 
it here as an illustrative check on our power-counting 
technique. 



B. The quantum phase ordering problem 

For the phase ordering problem, the right-hand side of 
Eq. (|2.5p gets multiplied by /(1/L,to ), 



dL 
~dT 



cA cOt 



mo). 



(3.3) 



In d = 3, the domain growth then displays four different 
power laws in different time or length regimes, as follows: 



L(t) oc < 



'(cA/i*)V3 t l/3 

(cQ L /L*)i 
(cALa/L*) 1 / 4 * 1 / 4 

-!/8 +1/4 
OC TOq t 1 



if L<{L l: L 2 } 
if Li < L < L 2 

if L 2 < L < Li 



(3.4) 



(cO L L 2 /L*)V2tV2 

octoq 74 ^/ 2 if {Li, L 2 }<L 



Compared to Eq. (|2.6p . the asymptotic time dependence 
of i remains unchanged, but the dependence of the pre- 
(3.1) factor on the equilibrium magnetization is m^ 4 instead 



of, 



1/2 



In the initial scaling regime, where L < L\, the 



time dependence is t 1 ^ 3 instead of W 4 in the classical 
case. In addition, there is an intermediate regime where 
L(t) grows as t if L\ < L 2 , and as t 1 / 4 if L 2 < L\. 

Equation (|3.4[) is the central new result of the present 
paper. It has been derived entirely by power counting, 
that is, from Eq. (|3.3[) which associated all gradients in 
the problem with powers of 1/L. Next we establish the 
validity of this procedure by means of a renormalization- 
group analysis that generalizes Bray's analysis of Model 



C. Renormalization-group considerations 

Adapting the renormalization procedure of Map we 
assign a scale dimension [L] = —1 to L, and a scale di- 
mension [t] = —z to time. Equation (|f ,ip suggests to 
choose the field 4>(x, t) to be dimensionless, [4>(x, t)] = 0. 
Let the Fourier transform of C in Eq. (|f ,ip be S(k, t) = 
J dx exp(—ik ■ x) C(r,t). The exponent r\ is defined by 
the structure factor S(k) = S(k,t — > 00) to behave as 
S(k) oc \k\~ 2+v ; this implies 77 = 2 — d. Position, time, 
and fields are then rescaled in the RG process accord- 
ing to x' = x/b, t' = t/b z , and <j)'{x',t') = <f>(x,t), re- 
spectively, with b the RG length rescaling factor. The 
free energy H, which has a naive scale dimension equal 
to zero, is assigned an anomalous scale dimension — y, 
H'[4>'] = b~ y H[<fi]. Finally, one needs to keep in mind 
that the functional derivative of H in Eq. (|2.ip removes a 
spatial integral and therefore acts, for scaling purposes, 
like an inverse volume with a scale dimension of d. A 
zero-loop renormalization of Eq. (|2.f p then yields 



Wn ,,W , ± , SH' 



(3.5a) 



with renormalized quantities 

A' = \b z - d - 2+ \ 

y = 1 b z - d+y . 



(3.5b) 



For Model B (7 = 0), the assumption that the trans- 
port coefficient A is not singularly renormalized at the 
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fixed point we are looking for (which assumes that a 
hydrodynamic description remains valid in the ordered 
phase), leads to the relation z = d + 2 — y,— which ex- 
presses the dynamical exponent z in terms of the energy 
exponent y. For Model J this fixed point is not stable: 7 
is relevant with respect to it. Assuming that 7 is not sin- 
gularly renormalized (which assumes that the spin waves 
in the ordered phase are characterized by w(fc) oc fe 2 ), 
leads to 

z = d — y. (3-6) 

The remaining question is the value of the anomalous 
energy dimension y. If defects in the order parameter 
texture determine the scaling properties of the energy, 
then y = d— 2 for a vector order parameter. 24 This yields 
z = 2, in agreement with the long-time behavior of L(t). 
For Model B, the same value of y yields z = 4= 4 With 
increasing length scale, one thus expects a crossover from 
z = 4 to z = 2, as is reflected in Eq. (|3.4[) . If L < L 2 , then 
one effectively has V rf ~ 2 in the free energy instead of V 2 , 
so one expects y = d — (d — 2) = 2. This leads to z = d 
for Model B, and z = d — 2 for Model J, as reflected in 
the first two lines in Eq. (|3.4p . The above considerations 
show that the naive power-counting considerations that 
lead to Eq. (|2.5[) or (|3.3p , which replace all gradients 
in the dynamical equation by 1/L, are indeed correct, 
subject to the above assumptions. Note that for an Ising 
order parameter, y — d — 1(21 and hence z = 3 for Model 
B, so the naive power counting breaks down. 9 

IV. DISCUSSION AND CONCLUSION 

In the remainder of the paper we provide a semi- 
quantitative discussion of Eq. (|3.4p by identifying the 
various length scales that enter the quantum phase or- 
dering problem. L\ has been identified in the context of 
Eq. (|2.6[) . L 2 can be identified from the explicit treat- 
ment of the ferromagnetic phase in Ref. |25J. We find 
L-2 = y/D/A, where D is the charge diffusion constant 
and A is the Stoner g ap or exchange splitting. A re- 
lated scale is Lt — y/D/T, which denotes the length 
scale where a nonzero temperature cuts off the nonana- 
lyticity. Finally, L* was introduced in connection with 
Eq. (|3.2p . Reference [l5| yields explicit expressions for 
the coefficients c and c, which give L* = ir£/72, with I 
the elastic electronic mean-free path due to the quenched 
disorder. Thus, 

Li = y/\/Sl h > L 2 = ^D/A, 

L T = VD/T , L* = ttI/72 . (4.1) 

We now estimate the values of these scales. J7l is on 
the order of emo/2m c c, with e and m e the electron charge 
and mass, respectively, and c the speed of light (not to be 
confused with the coefficient of the square gradient term 
in the Hamiltonian that is denoted by c everywhere else 
in this paper). The effects we are considering are largest 



if the magnetization is small; either because the system 
is a weak magnet, or because the quench is to just within 
the ordered phase (but outside the critical region; we 
consider a critical quench below). For a magnetization 
too = 10 - 100 G we have hn L /k B « 0.001 - 0.01 K. For 
the Stoner gap one expects HA/ks > T c , with T c the criti- 
cal temperature that corresponds to the parameter values 
after the quench. For low-T c magnets like MnSi or UGe2, 
this means HA/ks ~ 10 2 K. With free-electron parame- 
ters, and a Fermi wave number kp ps 1 A^ 1 , the mean- free 
path is related to the resistivity p by I ps 10 3 (/^ilcm/ p) A, 
and in the ordered phase one expects A ~ D ps hkp£/3m c . 
For p = 10 /iSlcm and T — 1 K, a rough estimate for the 
hierarchy of length scales thus is 

L* w 5 A , L 2 « 10 2 "A , 

L T w 10 3 A , L x » 10 5 A . (4.2) 

For L* < L < L 2 one has the initial i 1 / 3 behavior in Eq. 
(|3.4p . For L 2 < L < Li, the domain size will grow as 

L(t) cc to 1 ^ 8 i 1 / 4 , and for L > L\ the behavior crosses 
over to L(t) oc m^ 4 i 1 / 2 . With respect to the latter, one 
should keep in mind that domains larger than a few tens 
of microns are hard to achieve in zero magnetic field, 
except close to the critical points These predictions 
should be observable by time-resolved neutron scatter- 
ing. In particular, the magnetization dependence of the 
prefactor of the asymptotic t 1 / 2 law can be checked by 
quenching along trajectory A in Fig. [T] to different final 
pressure values. By recalling that the parameter c in Eq. 
(|2.2ap represents the square of a microscopic length scale 
that is on the order of an A we can estimate the time 
required for a domain to grow to sizes corresponding to 
the various length scales given in Eq. (|4.2p : 

t* fa 10~ 15 s , t 2 w 10~ n s , 

tr « 10~ 7 s , t x ps Is . (4.3) 

Notice that the microscopic time scale for the problem 
is given by the Fermi wave length divided by the Fermi 
velocity, which is about 10~ 15 s with free-electron param- 
eters. This is consistent with the value of t*. 

Now consider a quench into the critical region, which 
is divided into regimes denoted by I, II, and III in Fig. 
([1]). The classical critical fixed point controls Region 
I, where phase ordering has been discussed by Das and 
Rao! 3 - Regions II and III are controlled by the quantum 
critical fixed point, and the quantum critical behavior 
is known exactly^ Consider a quench to the quantum 
critical point, trajectory B in the figurei 2 ! The quantum 
ferromagnetic critical behavior is characterized by loga- 
rithmic corrections to scaling, which can be expressed in 
terms of scale dependent critical exponents. The dynam- 
ical critical exponent in d = 3 i a 6 ' 17 

z = 3 + const, x (In In b) 2 / In b, (4.4) 

to leading logarithmic accuracy. At the quantum crit- 
ical point, and in the context of domain growth, the 
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renormalization-group scale factor b represents 1/L(t). 
This leads to a growth law, with the time t measured in 
arbitrary units, 

L(t -> OO) oct l/3 e con a t.x(lnlnt) 2 _ (45) 

If the quench ends at a low temperature in the critical 
region, but not at the quantum critical point, L(t) will 
grow according to Eq. (|4.5|) until it becomes comparable 
to the correlation length £. In region lib, at longer times 
it will saturate at a value comparable to £, while in region 
Ha there will be a crossover to the asymptotic behavior 
as described by Eq. (|3.4p . A more complete description 
of critical quenches will be given elsewhere.— 

In clean systems, analogous mode-mode coupling ef- 
fects lead to a weaker nonanalytic term than in Eq. 
p.l[) ; in d = 3 it is fc 2 ln[fc|. However, the term is 
negative, which leads to a first-order transition^ The 
order of the transition is of no consequence for the 
phase ordering kinetics, but the requirement of a pos- 
itive transverse magnetic susceptibility in the ordered 
phase prevents the magnetization from ever being small 
enough for the nonanalytic term to dominate over the 
analytic one. For the magnon dispersion relation, one 
finds D{nio) = 7cmo[l — const, x ln(l/mo)], and the 
equation of state will ensure that D(mo) > 0j2& Sim- 
ilarly, for the phase ordering problem one has L(t) oc 
[1 — const, x \nt] t 1 / 3, in a transient regime, and L(t — > 



oo) oc mjj/ 2 [1 — const, x ln(l/mo)] t 1 ^ 2 asymptotically 
(const. > 0). The former result follows since, in the clean 
limit at T = 0, A(fe -> 0) oc l/\k\. 

We conclude by summarizing the original results ob- 
tained in this paper. First, we generalized the phe- 
nomenology that was developed to describe phase order- 
ing following a quench across classical phase transitions 
to the quantum phase transition case. Second, for the 
continuous quantum phase transition expected in disor- 
dered Heisenberg quantum fcrromagnets we obtained the 
growth laws in various time windows for quenches both 
deep into the ordered phase and to a point at or very near 
quantum criticality. Third, we gave the domain growth 
laws for clean itinerant Heisenberg quantum magnets, 
where the quantum phase transition is expected to be 
discontinuous. In the latter case the quantum effects are 
subleading due to the lower bound on the magnetization 
imposed by the first-order nature of the transition. All of 
these results are amenable to experimental verification. 
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